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ABSTRACT 

Non-uniqueness of deprojecting the integrated light distribution of a nearby or 
faraway triaxial body is reviewed in the context of deriving the volume density of the 
Galactic bar from the COBE/DIRBE maps of the Galactic plane. The exact origin of 
this non-uniqueness is studied. One can write down a sequence of triaxial bar models 
which appear identical in integrated light from the Sun's perspective, and the whole 
sequence is mapped out as a function of the Galactocentric distance of the observer 
from galactic to extragalactic distance scales. While mirror symmetries and perspective 
effects are compatible with any orientation of the bar in the sequence, weak upper and 
lower bounds can still be placed on the angles and axis ratios of the bar by positivity 
and other general requirements. Star count data of bulge giants are ideal for selecting 
a unique model from the sequence. 

Key words: Galaxy: centre - Galaxy: structure - Gravitational lensing - galaxies: 
photometry - galaxies: kinematics and dynamics 



1 INTRODUCTION 

Deprojection of galaxies from the observed light distribu- 
, tion on the sky plane to the intrinsic 3D volume luminosity 
distribution is one of the basic problems of astronomy. It is 
common knowledge that the deprojected results are gener- 
ally non-unique because of the freedom of distributing stars 
along any line of sight. General constraints such as sym- 
metry, positivity and self-consistency are often too weak to 
completely break the degeneracy of models, in particular, 
the fundamental degeneracy of a triaxial ellipsoidal distribu- 
tion vs an axisymmetric spheroidal one; both can project to 
the same elliptical isophotes in the line-of-sight integrated 
light distribution (Contopoulos 1956). Kinematic informa- 
tion such as rotation around the apparent minor/major axes 
of the isophotes and the 2D projected velocity field is nec- 
essary to solve this so-called "shape problem" of elliptical 
galaxies (Binney 1985, Franx, Illingworth & de Zeeuw 1991, 
Statler 1994 and references therein). 

Nevertheless, as Blitz & Spergel (1991) pointed out in 
the case of the Galactic bar vs. an oblate bulge, the degener- 
acy is at least partially broken by photometric information 
from the perspective effect of nearby triaxial or bi-symmetric 
objects, namely with mirror symmetry with respect to three 
orthogonal planes. Perspective is the basis for the detection 
of the Galactic bar in the maps of the Galactic plane made 
by the DIRBE experiment aboard the COBE satellite. The 
dust-corrected infrared maps from COBE show significant 
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systematic brightness variations between the right hand side 
(/ > 0°) and the left hand side (I < 0°), and small system- 
atic differences between the lower (b < 0°) and the upper 
(b > 0°) Galactic plane (Weiland et al. 1994): a strong signa- 
ture of a nearly edge-on (mirror symmetry with the b — 0° 
plane) triaxial bar with the major axis of the bar rotated 
by an angle from our line of sight. A simple sketch of the 
geometry is shown in the top left diagram of Fig.l). 

Dwek et al. (1995) made the first attempt to fit the 
COBE/DIRBE maps with a few simple analytical luminos- 
ity models. They found that the COBE maps are grossly 
consistent with several triaxial models: a very strong bar 
pointing almost along the Sun-Galactic center line of sight 
or a short bar pointing more sideways. Despite still signifi- 
cant residues of all their models, the existence of reasonable 
fits of the COBE map with a wide range of bar angles, axis 
ratios, radial profiles and degrees of boxyness suggest an in- 
trinsic non-uniqueness among triaxial models which is not 
lifted by the perspective effect. 

Binney & Gerhard (1996) separated the non-unique ori- 
entation of the bar from the likely unique 3D distribution 
once the three mirror planes of the bar are fixed. Their nu- 
merical experiments show that the final result of their non- 
parametric deprojection depends only on the assumed ori- 
entation of the bar, but not on the initial seed models for 
starting their Lucy-Richardson iterations. They explain the 
latter uniqueness as a compromise between the symmetry 
requirements of the density and a good fit to the perspective- 
induced asymmetries in the COBE map. 

Binney, Gerhard & Spergel (1997) applied this non- 
parametric method to the COBE map after correcting for 
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dust with an advanced extinction model from Spergel, Mal- 
hotra & Blitz (1997). They concluded that the major axis 
angle of the bar is likely very close to 20° with a good fit to 
the COBE map still possible for a positive smooth triaxial 
bar with major axis angle between 15° and 35°. They sus- 
pect that these limits on the bar angle mark the boundary of 
unphysical densities with negative regions and/or peculiar 
profiles. 

Compared to external systems, perspective provides the 
unique signature for the Galactic bar as well as the subtle- 
ness to its non-uniqueness. The non-uniqueness in a nearby 
bar is entangled with perspective effects, making it a to- 
be-defined narrow sub-class of non-uniqueness for external 
systems. Simple transformations such as shearing/stretching 
an external ellipsoid in the line-of-sight direction can gener- 
ate a sequence of ellipsoids with the same surface brightness 
map (Stark 1977, Franx 1988). But as Binney & Gerhard 
(1996) emphasize, moving material along the line of sight 
is forbidden by the mirror symmetries of a nearby bar. The 
well-known trick with axisymmetrical galaxies of adding a 
small amount of unphysical density with zero surface bright- 
ness (Rybicki 1986, Palmer 1994, Kochanek & Rybicki 1996, 
van den Bosch 1997, Romanowsky & Kochanek 1997), chris- 
tened konuses by Gerhard & Binney (1996) for extragalactic 
axisymmetric bodies would not apply either. While preserv- 
ing the orientation and the symmetry of the bar, adding 
konuses perturbs the left-to-right difference map. 

Unfortunately none of the previous analyses with the 
COBE bar and external bars offer a satisfactory explana- 
tion for the non-uniqueness. In particular, the true source 
of non-uniqueness is never clearly singled out, and how this 
non- uniqueness varies as a function of the distance to the 
body is never clearly stated. It is also unclear how to trans- 
form one solution to another and map out the sequence of 
models admitted by the COBE maps, which is one of the 
main sources of information on the stellar distribution of 
the Galactic bar. Answers to such questions are important 
for integrating the constraints from the COBE map to other 
data of the bar from star counts, gas/stellar kinematics, and 
microlensing (see de Zeeuw 1993, Gerhard 1995). For exam- 
ple, the angle of the bar is a main source of uncertainty 
in comparing the COBE map with observed microlensing 
optical depth of the bar (Zhao & Mao 1996, Bissantz, En- 
glmaier, Binney & Gerhard 1997); it is also a main source 
of uncertainty in determining the number density of brown 
dwarfs in the Galactic bar using microlensing (Zhao, Rich 
and Spergel 1996). 

This paper gives an analytical description of the nature 
of the non-uniqueness in the nearby COBE bar and points 
out its main differences with non-uniqueness in external sys- 
tems. I also show how the positivity of the distribution puts 
a loose limit on the bar models, and discuss the best meth- 
ods to lift the remaining degeneracy. 



2 THE NATURE OF NON-UNIQUENESS 

In this section I will describe the sources of non-uniqueness 
purely in words and illustrations without getting into the 
more mathematical aspect of the problem, which will be 
supplimented in §3 together with some applications to the 
COBE bar. 



2.1 Mirror symmetries vs. invisible densities 

Consider for the time being only bars which are symmet- 
ric with respect to the 6 = 0° plane. Divide the Galactic 
bar and the COBE map into a left and a right part with 
the I = 0° plane, which passes the Sun-center line and the 
rotation axis of the Galaxy. When folded along the / = 0° 
line, the COBE/DIRBE Galactic plane map 1(1, b) can be 
decomposed into two independent maps: a difference map 
1(1, b) — I(—l,b) by subtracting the I < 0° side from the 
I > 0° side, and an addition map 1(1, b) +I(—l,b) by adding 
up the two sides. These two maps are line-of-sight integra- 
tions of the odd and the even parts of the volume density 
v of the inner Galaxy, v dd and v even , decomposed in terms 
of the symmetry with respect to the I — 0° plane. The dif- 
ference map and v dd are anti-symmetric with I = 0° plane, 
the addition map and ^ eue „ are left-to-right symmetric. 

Non-uniqueness originates from the large degrees of 
freedom in fitting the left-to-right symmetric addition map 
1(1, b) + I(— I, b) |j any even component, such as an axisym- 
metric bulge or disc, a long end-on bar (the major axis co- 
incides with the Sun-center line), a short side-on bar (the 
middle axis coincides with the Sun-center line) look very 
much alike in projection. In fact, in subtracting the axisym- 
metric component from the end-on bar model, one makes an 
unphysical distribution with some negative density regions, 
which is similar to extragalactic konuses in the sense that 
its projected light intensity is exactly zero after integration 
over any line of sight. I shall call such left-to-right symmetric 
unphysical components "invisible densities" since superim- 
posing them on the bar is undetectable in projected light. 
They are the sources of non-uniqueness. 

An example of an invisible density is illustrated in the 
upper right corner of Fig.l, where I tailor the radial profile 
of a stratified ball-shaped bulge in such a way so to have the 
same angular size and projected intensity as a uniform cigar- 
shaped bar placed end-on. Subtracting the ball-shaped bulge 
from the cigar-shaped bar yields an invisible density. The 
thin dark ring in the bulge is a density peak, corresponding 
to the line of sight to the far edge of the cigar-shaped bar, 
a direction where the depth and the projected intensity of 
the bar are at maximum. The fall-off of density towards the 
edge of the bulge corresponds to the ever-decreasing depth 
of the cigar-shaped bar with increasing impact parameter of 
the line of sight. Now superimpose a finite fraction, say /, 
of this unphysical component on the COBE bar, say, it is an 
ellipsoidal bar with a Gaussian radial profile. While such an 
"operation" is invisible in the projected light, its result is a 
new density with generally so irregular isophotes that there 
is no specific symmetry plane. 

When the mirror symmetries with respect to three given 
orthogonal planes are imposed on the density distribution 
of the bar, the forementioned "invisible operation" is nearly 
forbidden because of the general twist of the major axis po- 
sition, but not always. The exception is, as shown in the 
enlarged lower diagram of Fig.l, when the superimposed 
end-on bar is a clone of the original bar, in which case the 



t Deprojection of the difference map 1(1, b) — I(—l, b), which dis- 
tinguishs a bar from an axisymmetric bulge, is assumed to give a 
unique odd solution u dd in this paper. 
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Figure 1. Crossections of a uniform cigar-shaped bar tilted at an 
angle a from the Sun-center line (upper left), a two-component in- 
visible density (upper right), and a superposition of the above two 
bodies (the enlarged lower diagram); several line-of-sight pathos 
from the Sun are drawn in thin solid lines. The invisible density is 
made by subtracting from an end-on cigar-shaped bar a stratified 
ball-shaped bulge of the same angular size and projected inten- 
sity. Adding such an invisible density to the cigar-shaped bar on 
the upper left has no effect on the latter's projected intensity, but 
rotates the mirror plane of the spatial distribution to the dotted 
line in the middle at an angle ^ from the Sun-center line. 



final density should have mirror symmetry with respect to 
a new plane (the dotted line) which divides the angle of 
the Sun-center line and the major axis of the original bar 
in half. Namely the mirror symmetry is preserved, only the 
symmetry plane is different. Subtracting the round bulge 
has no further effect on the symmetry plane, but will take 
away any trace of transformation in the addition map of 
the integrated light. So an invisible density which rotates 
and preserves mirror symmetry is obtained by placing the 
original bar at end-on then subtracting off a round bulge. 
Generally the new major axis position depends on the frac- 
tion / of the added invisible density. When / is sufficiently 
small, the new density could also be positive everywhere. 
The general technique of designing invisible densities for a 
realistic triaxial bar is shown in §3.1 and §3.3, in particular, 
how to rotate the bar without twisting its major axis with 
radius. 



2.2 Perspective effect, extinction and best-fit 
models 

How is the perspective effect of the bar decoupled from 
the non-uniqueness problem? The perspective effect is only 
prominent in the difference map, and is by-passed by re- 
stricting all the transformations to the even part of the den- 
sity: v dd is never affected by the invisible density. Still the 
invisible density is a function of the Galactocentric distance 
of the observer Ro- For a hypothetical observer receding 
from the Galactic bar towards the anti-center, its Galac- 
tocentric distance Ro increases, (the axisymmetric part of) 
the invisible density should be adjusted slightly to accomo- 
date the changing perspective so to preserve this kind of 
non- uniqueness all the way to extragalactic distance. The 
common ambiguity of an extragalactic axisymmetric bulge 
with an end-on bar is a very special case of the kind of non- 
uniqueness with the COBE bar. As 1 will come back to in 
§3.2, other distinctive features of the invisible density are 
that it has a net luminosity, it is not confined to any cone 
in the Fourier k-space and it is symmetric with respect to 
the I = 0° plane. 

Will invisible density remain invisible in the dusty 
Galaxy? Even at the near infrared the extinction of the 
COBE bar is a strong source of uncertainty for deprojection 
as shown by the strong color variation in the COBE/DIRBE 
maps (Arendt et al. 1994). In the first order approximation 
where dust is primarily situated in the foreground within the 
near 2 kpc of the line of sight to the center (Arp 1965), the 
projected intensity of the bar is simply multiplied by a red- 
dening factor, thus the zero projected intensity of the invisi- 
ble density is not affected. However at low latitude (|6| < 3°) 
where dust and stars are mixxed both in the disc and in the 
bar regions (Spergel, Maholtra & Blitz 1997), extinction will 
modulate the bar density with an ill-determined function for 
the dust distribution, which should be deprojected simulta- 
neously with the stellar distribution. The invisible density, 
defined in the dust-free context, will cast faint (irregular) 
shadows on the COBE/DIRBE maps; one sees the near end 
of the end-on cigar-shaped bar (cf. upper right of Fig. 1) 
faintly in a heavy mist of interstellar dust. Although in prin- 
ciple one could still define a certain density which projects 
to zero after being modulated by dust, and infer the dust 
distribution from, say, star count data, independent of the 
to-be-deprojected maps, but finding such a density which 
preserves the triaxial symmetry of the bar becomes hope- 
lessly difficult. 

Do "best-fit" bar models mean anything? Since for any 
bar which fits the COBE/DIRBE maps, one can build a se- 
quence of triaxial models with any bar angle to fit the data 
with exactly the same accuracy by adding different amount 
of the invisible density, it is ill-dehned to speak of the "best- 
fit" bar model or bar angle from the COBE/DfRBE maps. 
The "best-fit" bar angles derived in the parametrized models 
of Dwek et al. are model dependent determinations. These, 
which could have fairly small error bars and/or residues (e.g. 
their Gl, G2 and E3 models), are a mere reflection of the 
different ranges where their different a priori assumptions of 
the boxiness and radial profile of the density are "best" satis- 
fied. The systematic residues in their models reflect only the 
unphysicalness of the a priori assumptions in detail (their 
diamond shaped bar is the best example). In fact it would 
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not even be surprising to have a perfect fit to the data for ar- 
bitrary bar angle (cf. eqs. |22[ ^) provided the dereddened 
COBE maps are sufficiently smooth because there are in 
principle enough degrees of freedom in the mathematical 
functional forms of a triaxial density to make the problem 
underconstrained. It is straightforward to show that for ob- 
servers nearby or at infinity this kind of ambiguity with the 
major axis angle occurs in every even m mode (m = 2, 4, ...), 
odd m mode (m = 1,3,...), and m-fold spiral arm mode 
where the major axis angle a also changes with R, and also 
in any band-limited superposition of finite number of modes. 
In all cases the major axis angle cannot be determined solely 
on the basis of the goodness of fit to the COBE map. 

Then where do the limits for the bar angle come from? 
The only hope left to constrain the bar angle from the COBE 
maps is to impose some generally sensible unbiased require- 
ments of the density, such as positivity, and to some extent 
smoothness, regular shape and radial profiles, as in Binney, 
Gerhard & Spergel (1997). An example is given in the fol- 
lowing to illustrate the main points in this section. 



3 APPLICATION TO THE COBE BAR 

So far I have assumed that one can design invisible densities 
which project strictly to zero even for nearby bodies, and I 
have also restricted to the simple case that the observer is 
in the mid-plane of the bar. In this section some examples 
of the invisible density are written down as explicit func- 
tions of the observer's distance and viewing angles to the 
bar, allowing for a small vertical offset of the Sun from the 
mid-plane of the bar (Dwek et al. 1995, Binney, Gerhard 
& Spergel 1997). These toy models are used both to illus- 
trate the non-uniqueness of the COBE bar and and to study 
possible methods to break the degeneracy. 



3.1 An m — 2 bar viewed from a general distance 
and perspective 

An m = 2 triaxial density distribution with any orienta- 
tion of symmetry planes can be fairly generally written in a 
Galactocentric polar coordinate (r, 8, <fi) as 



f/=o(r,0,< 



G(r) + P(r)H f=0 (e, 



(1) 



which is a superposition of a spherical component G(r), 
and a triaxial perturbation with the radial profile P(r) and 
the angular part Hf = o(8,<j)); the latter equals to a function 
Hf(0, (j>) evaluated at / = 0, where 

2 n 



H 



s + (/ + Co) cos 
+2(ci cos 4> + si sin <p) sin 8 cos ( 
+ (c2 cos 2(f> + S2 sin 2<j>) sin 2 8. 



(2) 



One can show that the perturbation is triaxial since 
r 2 Hf=o(8,<f>) reduces to a quadratic function in rectangu- 
lar coordinates with ellipsoidal isosurfaces whose symmetry 
planes are determined by the constants Co, ci, C2, so, si, and 

S2- 

Now a sequence of triaxial models with the same pro- 
jected density as i/f = o(r, 8, (f>) can be made if one can write 
down its invisible density, which should be made by sub- 
tracting a spherical model from an end-on prolate bar with 



amplitude P(r). For convenience fix the pole of the Galac- 
tocentric polar coordinates (r, 8, <j)) as the Sun-center line of 
sight, and the zero point of the roll angle <j> as the 6 = 0° 
plane, so that in this coordinate system P(r) cos 2 8 would de- 
scribe an end-on prolate (i.e., (^-independent) distribution. 
Now suppose a spherical density S(r) can be inverted* from 
the integral equation, 



S(r)dD 



[P(r) cos 2 8] dD = T(l,b), 



(3) 



where J^°dD is an integration along any line-of-sight direc- 
tion (l,b), then 

vj (r, 8, 0) = v f=0 (r, 0, 0) + [P(r) cos 2 6 - S(r)] f, (4) 

describes a sequence (as function of the parameter /) of tri- 
axial models with the same projected density. The sequence 
clearlypreserves triaxial symmetry because like Vf=o(r, 8, (j>) 
of eq. [l| 

v f {r,8A) = [G(r) - fS(r)] + P(r)fl>(M), (5) 

also consists of a spherical component G(r) — fS(r), and a 
triaxial perturbation with an angular dependence Hf(8,(f>) 
similar to Hf=o(9,(j)) (cf. eq. ||) except that the orientation 
of the symmetry planes will also depend on /. 

By tuning solely / the model (cf. eq. |H|) changes perspec- 
tive with the Sun: the offset of the Sun from the mid-plane 
and the rotation of the major axis in the bar mid-plane 
are varied together along a ID sequence. The less interest- 
ing third angle which describes a roll of the bar mid-plane 
around the Sun-center line, and is equivalent of the sky posi- 
tion angle in the case of extragalactic bodies, is not affected 
by adding an invisible density; the latter is invariant of a 
roll around the line-of-sight axis. 

As a specific example to show how the model changes 
with the observer's distance and perspective, let 

r 2 r 2 r 2 

G(r) = v oex p(- — ), P{r) = v —exp{ — -), (6) 

where a is the characteristic scale of the model. The spherical 
density S(r), as specified in the integral eq. ^| can be solved 
analytically with a variation of the well-known Eddington 
formula for deprojecting spherical system. 0The dependence 
on the Galactocentric distance Ro is generally such that if 
Sao(r) is S(r) in the limit of _Ro — > oo then S(r) — Soo(r) oc 
(^■) • For the models in eq. g, 



S(r) = Soc(r) 



Rl 



2r 4 



3r< 



■exp( 



,(7) 



t Generally if J(p) is the integrated intensity of the light 
distribution S(r) or P(r) cos 2 8 along a path with a Galacto- 
centric impact parameter p, including the contributions from 
both the forward direction and the backward direction, T(l, b) 
and T(tt + 2,-6), then J(p) = T(l, b) + T(n + 2,-6) = 
cos 6du, where u is the offset along 

the line of sight from the tangent point, r = yjp" 2 + u 2 , and 

r£2sl = {^^Jl - j£ - . S(r) is then inverted with an Abel 

transformation S(r) = — 7r _1 J°°— jf^~{P 2 ~ r 2 )~ 1 / 2 dp. 
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Figure 2. Equatorial slices of two triaxial models with identical 
surface brightness map, one at a = 25° (thick dashed contours) 
and the other at a = 50° (thin solid contours). The Sun is at 
(—1, 0) to the central left. Two solid lines are drawn to show a 
line of sight at positive and negative longitude respectively. 



and 
v s (r, 



2 2 

= ^oexp(-^){l^e X p(-^)[ 



(8) 



2 +f m I ~ 2 i 



where the extra damping factor exp( — js) in front of the 
square bracket is to guarantee that the model is spherical 
and positive at large radius. The above equation yields tri- 
axial models with monotonic and positive density and vari- 
ous degrees of flattening/triaxiality and boxyness for a wide 
range of scale length a and shape and perspective parame- 
ters (co, Ci, C2, so, si, S2) and /. The models resemble the Gl 
and G2 models of Dwek et al. (1995) for the Galactic bar 
in terms of a Gaussian radial profile. Since the Sun is only 
about 10 pc from the midplane of the bar (Binney, Gerhard 
& Spergel 1997), the parameters si and S2 should be of the 
order 10pc/7?o ~ 10/8000 ~ 0.01, much smaller than unity. 
In particular, models with 



si = S2 — 0, / = C2 — cq — 2ci cot 2a, 



(9) 



prescribe a sequence of edge-on (mirror symmetry with b = 
0° plane) m = 2 perturbation, where a is the angle between 
the major axis of the perturbation and the I = 0° plane 
which passes the Sun and the rotation axis of the Galaxy. 
Fig. 2 shows an equatorial slice of two such edge-on models 
from the sequence at a = 25° and 50° respectively. The 
shape of bar is a function of the major axis angle such that 
models in a sequence yield the same projected light intensity 
from the Sun's perspective. 



3.2 Non-uniqueness of the COBE bar vs. external 
systems 

Now I digress from the COBE models to stress a few key 
differences between the non-uniqueness of the COBE models 



and external systems. Unlike its external counterpart with 
Ro — * 00, the invisible density for a nearby triaxial body is 

not massless. For the bar models here, the invisible density 

2 

has a mass proportional to The total luminosity of the 



bar models (cf. eqs. ^ and ^) 



d rv(r, 



TV * VQCL 



2' + 



3so + c 3 a 
2 2 ~Rj 



,(10) 



changes with /. 

To fit the same projected light intensity the models 

should have the same normalization vq, but the total lu- 

2 

minosity changes by a factor on the order of ^ ~ a few 



percent for the COBE bars. This is because unlike external 
systems, the observed angular distribution of light in the 
COBE bar does not simply sum up to a unique measurement 
of its total luminosity. The result depends on knowing the 
distance to the bar and the orientation of the bar. The upper 
panel of Fig. 3 shows that an observer at the Sun's position 
i?o = 8kpc may slightly underestimate/overestimate the in- 
trinsic luminosity by about 10% depending on assumptions 
of the bar angle. Moving out to the anti-center the observer 
has a full outside view of the bar and hence a better deter- 
mination of its luminosity but sacrifices somewhat tighter 
constraint on the bar angle from the perspective effect and 
positivity. Going to the center of the bar, the observer has 
full information of the orientation of the bar (the brightest 
direction is the major axis direction), but very poor informa- 
tion of the intrinsic luminosity profile and total luminosity. 

Another feature of the invisible density here is that un- 
like konus-like structures, it is not confined to any cone in 
the Fourier k-space even in the limit Rq — > 00. This is shown 
explicitly by computing the Fourier transform of the (pro- 
late band-limited) invisible density (cf. eqs. ^| and [?]) using 
eqs. 4, 24 and 25 of Palmer (1994). 



F(k) = J d 3 r exp(-ik ■ r) [P(r) cos 2 - S(r)] 



= F fl0 ^ oo (k) + - ff T(fc) ) 
tt 



(11) 



(12) 



which is a k-space prolate distribution around the Sun- 
center line, where 



Fr 



3 3 / 2 2\ 
^oo(k) = 7r2 v a Icos k — -J 



2\ k 2 a 2 k 2 a 2 
cos k - -J — p exp( — ), (13) 



0k is the angle of the k-vector with the Sun-center line, and 

m ,, , 2 3/3 21fcV fc 4 a 4 \ . fcV, 
T(k) = 7r 3 Voa \ -- + — —J exp(-— ). (14) 

One recovers the mass of the invisible density in the limit 
k -» 0, F(0) = f^T(0) = -f^OrtiW) (cf. eq[fij. Eq. g 

shows that even when the bar is placed at infinity, F(k) = 
Jfl ->°o(k) 7^ everywhere (except at k= and/or 0k = 
cos -1 v/^). Since F(k) = for konus-like structures (or 
its triaxial version) outside a certain "cone of ignorance" 
around a principal axis of the model (Rybicki 1986, Gerhard 
& Binney 1996, Kochanek & Rybicki 1996), the kind of non- 
uniqueness shown in this paper has little to do with konuses. 
The two sequences of models meet only for an end-on prolate 
bar seen at infinity, which can alternatively be interpreted as 
an axisymmetric system with inclination i = (face-on). In 



© 1997 RAS, MNRAS 000, §-| 



6 



this special case the invisible density here becomes a konus 
but the half opening angle of the "cone of ignorance" is as 
large as ^. 

The invisible density is also of different nature from 
the kind of non-uniqueness associated with a simple shear 
and/or stretch transformation of an ellipsoid. This applies 
even in the limit that the object is at infinity because these 
transformations normally change the odd part of the density, 
while the invisible density here is strictly an even function. 



3.3 Limits on the bar angle vs. positivity 

The limits for the bar angles are given by the positivity re- 
quirement. They can be derived analytically in the context 
of models more general than eq. ^ without additional effort. 
As it turns out, the results are simplest for the broad class 
of interesting edge-on models with an axisymmetric com- 
ponent and an m-mode perturbation; m = 2,4, ... if mirror 
symmetry is imposed. These models have a density distri- 
bution which can be written with no loss of generality as 



v a (R,Z,ip) = 

A a (R,z) 
K m (R,z,ip) 



A a (R, z) + 5m ^' Z * cosm(V> - a) 



C(R, z) — cot maV m {R, z), 
[B m (R, z) cos mtp — V m (R, z)\ . 



and 



dDK m (R,z,iP) = 0, 



(15) 

(16) 
(17) 

(18) 



where a Galactocentric cylindrical coordinate system 
(R, z, ip) is used with the azimuthal angular coordinate ip 
counted from the I — 0° plane; the major axis of the per- 
turbation is along ip = a with amplitudes Emilhll w here 

° ' r sin ma ' ■ . 

B m (R,z) specifies the odd part of the density (cf. eq. |23| ), 
and C(R, z) is an axisymmetric fit to the addition map (cf. 
eq. ^), and K m (R, z,ip) is an invisible density, made by 
subtracting an axisymmetric component V m (R, z) from an 
end-on bar B m (R, z) cos mip. 

A few brief comments before proceeding to derive the 
limits on the bar angle. Eqs. jl5|- [l^ imply that any triaxial 
bar with the same amplitudes C(R,z) and B m (R,z) will 
follow a sequence with identical line-of-sight integrated light 
distributions on both positive and negative I sides I(±l,b), 
the same odd parts v dd (with respect to the ip = / = 0° 
plane) in density, and different even parts u everl , but with the 
difference being a linear superposition of invisible densities 
K m (R,z,ip). For triaxial densities with generally multiple 
modes (m = 2,4, ...), they are related to each other by an 
invisible density 



E 



(cot ma — cot ma)K m (R, z, ip), 



(19) 



where the sum of different modes in the given proportion 
is necessary to prevent spiral-arm-like twist of the major 
axis. Adding this invisible density to a model with symmetry 
plane at ip = a rotates the plane to ip — a' . The sequence 
of edge-on models given by eq. [| and ^ is a subclass of the 
models here with m = 2, and 

R 2 

B 2 {R,z) = - Cl P(r)— , V 2 (R,z) = -2 Cl S{r)-B 2 (R,z),(20) 



C(R,z) = G(r)+(c -c 2 )S(r)+P(r) s - c 2 + 2c 2 



R 2 



(21) 



Generally C(R, z) and B m (R, z) should be derived from fit- 
ting the addition and the difference maps of the projected 
light, 



[I(l,b)+I(-l,b)] = 



[!(/,&) -!(-/,&)] 



C(R,z)dD, 



(22) 



B m (R,z)sinmipdD, (23) 



and V m (R, z) from eq. |l7j . 

Impose positivity requirements to v a (R,z,ip), 



A a (R, z) = [C(R, z) - cot maV m (R, z)] > 



B m (R, 



(24) 



for all (R,z). This reduces to an upper and a lower bound 
for the angle a, 



Max [t_] < tan — < Min [t+] , 



(25) 



where t_ < t < t + is the range bounded by the effectively 
quadratic inequality for a variable t 



B„ 



Vrr, 



B m + Vrr, 



2C L *+ "'" 2C "' t * <1 (26) 

at a given position (R, z) on the meridional plane, and the 
overlapped interval of these ranges is used in eq. ^5| The 
exact upper and lower bounds of the angle a for the COBE 
bar will be discussed elsewhere as it involves solving three 
integral equations numerically (eqs. ^ and ^ for B m (R, z) 
and C m (R, z) and eq. [l?] for V m (R, z)), and several practical 
issues not considered so far (including a detailed dust model 
for the COBE/DIRBE maps, the choice of fitting regions 
with reliable photometry, the degree of smoothness of the 
solution in a non-parametric fit, and search of solutions in 
the ID parameter space spanned by a). Here it is only worth 
commenting on two points. First Binney, Gerhard & Spergel 
(1997) showed that it is possible to obtain good fits to the 
COBE maps with a triaxial bar if the bar angle is in the 
range 15° < a < 35° and the Sun is very close the midplane 
of the bar. A sequence of models which is roughly consistent 
with these limits is discussed in the next section. Second 
the range for a is set by C m and B m , which are in turn 
set by the addition map and the difference map resp ectively 
with an integration over the line of sight (cf. eq. |i^and ^) . 
Hence the range for a is also a function of the distance 
to the bar. However, the upper and lower limits are only 
weakly dependent on the distance as long as the observer is 
outside the bar, as shown in Fig. 3 by the nearly horizontal 
boundaries of the shaded region for Ro 3> 2 kpc. 



3.4 A toy model and implications to the COBE 
bar 

To illustrate the non-uniqueness in the Galactic bar, and 
how to lift the degeneracy, consider edge-on models in eq. ^ 
and ^| with 

a = 2 kpc, c = c 2 = -y = i, so = -1, (27) 



and (cf. eq. []) 



si = s 2 = 0, / = C2 — Co — 2ci cot 2a — — cot 2a, (28) 
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Figure 3. For an edge-on bar with a fixed M/L and parameters 
given in eq. E7I, positivity of the densi ty p uts loose upper and lower 
limits on the bar angle (cf. eqs. [25) , [26^ |2oj and ^j). These limits 
are functions of the distance to the bar, which mark the upper 
and lower boundaries of the shaded region in the log(i?o/8kpc) 
vs. a plane. The uncertainty of the bar angle propagates to the 
inferred value for the luminosity of the bar. The contours show the 
bar's absolute magnitude (upper panel) in the log(i?o/8kpc) vs. a 
parameter plane with an interval of 0.05 magnitude, and the equal 
contours of the central escape velocity -y/ 2 1 <J>(0) (lower panel) 
in intervals of 5%. Both the escape velocity and the luminosity 
increase with increasing a, but the increment is small as long 
as one is observing from well outside the bar, despite the strong 
ambiguity with the angle of the bar. 



which prescribes a sequence of edge-on m — 2 bars with a 
Gaussian radial profile inside about 4 kpc of the center with 
major axis at an angle a from the Sun-center line. The den- 
sity is strictly positive everywhere for models with a broad 
range of major axis angles, 22.5° < a < 63.5° for Ro = 8 kpc 
(cf. the shaded region in Fig. 3). Models are also monotonic 



everywhere except for a small bump near 



4 kpc on 



the short and middle axes. These limits will be significantly 
relaxed when an axisymmetric exponential disc is included 
because it adds a positive background distribution; a disc 
with peak intensity about 1/10 of the bar can also make 
a ~ 15° and a ~ 75° models positive everywhere. 

While the specific models here are clearly toy models, 
too simplistic to fit the COBE map in any detail, as shown 
in Fig. 2, they grossly resemble the Galactic bar in terms of 
its Gaussian radial profile (Dwek et al. 1995), its shape and 
axis ratios. In particular, the range of the bar angle here 
is largely in agreement with those in the literature (see a 
table compiled in Gerhard 1995), including the results of 



Binney, Gerhard & Spergel (1997) and Dwek et al. (1995) 
derived from fitting the COBE map, and Stanek et al. (1994, 
1997) derived from fitting the star count data of bulge red 
clump giants. Hence the models are useful to demonstrate 
the non- uniqueness and test various techniques for lifting 
the non- uniqueness. 

Fig. 2 shows an equatorial slice of the model at a = 25° 
(the thin peanut-shaped bar) in heavy dashed contours and 
a model at a = 50° (the more concentrated bar) in solid 
contours. These two grossly different models look exactly 
the same in line-of-sight integrated intensity for an observer 
at (—1,0) at the central left (the assumed position of the 
Sun). Neither of these two models can be easily ruled out 
on the basis of stability and self-consistency. 

Nevertheless it is possible to break the degeneracy with 
kinematic data since the invisible density changes the ori- 
entation of the bar potential model as well as its shape and 
depth. For the models in eq. |^, the depth of the potential 
well at the center is 



|$(0)| = irv a 2 G(M/L) 



2s + -co 



- — 

3 Rn 



f 



,(29) 



which decreases linearly with increasing fraction (/) of the 
superimposed invisible density, where M/L is the mass-to- 
light ratio, assumed to be the same for all models. As a re- 
sult, the potential well of the a — 25° model is slightly shal- 
lower than that of the more centrally concentrated a = 50° 
model (cf Fig. 2). However, the difference is only 7% in term 
of the maximum escape velocity v/2|$(0)| of the models 
(cf. lower panel of Fig. 3). The differences in terms of the 
mass weighted average velocity dispersion (estimated from 
the virial theorem) and the circular velocity (estimated di- 
rectly from the potential) are also at only a few percent 
level, too small to be measured with certainty. To distin- 
guish the models with kinematic data, one needs to model 
the anisotropy and perspective effects of the bar orbits in 
detail, e.g., the longitude- velocity diagram for orbits of gas 
clouds (Binney et al. 1991, Zhao, Rich & Spergel 1996), and 
the vertex deviation in the stellar velocity ellipsoid (Zhao, 
Spergel & Rich 1994); these models introduce additional 
complexity such as the pattern speed of the bar, and the 
stellar distribution function. 

The most straightforward approach to break the degen- 
eracy of the COBE bars is perhaps to compare the models 
with measurements of the line of sight distance distribution 
of the bar from color-magnitude diagrams. The depth of the 
models, defined here as the rms dispersion of distance in the 
Sun-center line of sight, decreases from 0.26i?o for a = 22.5° 
to 0.19i?o for a = 63.5°. The 30% difference is detectable 
with good distance indicators and large data sets, such as 
the bulge red clump giants in the microlensing and variable 
star surveys towards the Galactic bulge (Paczyriski 1996). 
It is also interesting that the invisible density also leaves 
observable traces on the microlensing maps, i.e., the distri- 
bution of the optical depth on the sky plane (Evans 1995, 
Bissantz, Englmaier, Binney & Gerhard 1997). This is be- 
cause the optical depth is proportional to both the projected 
intensity and the line-of-sight depth, and while an invisible 
density does not change the former, it changes the latter. 
The microlensing optical depth increases by about 50% by 
going from the a = 63.5° model to the a = 22.5° model. 
Nevertheless, since the data set on bulge red clump giants 
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Figure 4. Simulated distance modulus distributions for red 
clump giants in the two bar models a = 25° (dashed curves) and 
50° (solid curves) in Fig. 2 at line of sight directions (±5°, —4°). 



is on the order of a million times larger than that of mi- 
crolensed sources, it argues for the former plus the COBE 
map being the speedest approach. However, fitting the mi- 
crolensing optical depth, gas/stellar velocities can set the 
normalization of the density, namely the mass-to-light ratio. 
To fit the same projected light intensity the models should 
have the same normalization uq% but the total luminosity 
increases by about 8% (cf.eq. [10) by going from a a = 22.5° 
model to a a = 63.5° model. 

Fig.4 shows the simulated star count data for the two 
models at line of sight directions (I, b) = (±5°, —4°), with an 
assumed Galactocentric distance modulus (DM) 14.5 mag- 
nitude. The predictions have taken into account the growing 
volume with radius for a given solid angle, and have been 
convolved with the dispersion of the absolute magnitude of 
the clump giants (about 0.2 magnitude). The predicted dis- 
tance modulus distributions are quite different for the two 
models in terms of the peak position, the width and shape 
of profiles. The fact that for both models the I > side is 
closer than the / < 0° side is because the major axes of both 
bars are placed at positive a. Interestingly the amount of 
I > 0° vs. / < 0° asymmetry is the same for both models, as 
indicated by the two short vertical lines in Fig.4, because the 
left-to-right symmetric invisible density has no effect on the 
asymmetry. The distributions resemble Fig. 2 of Stanek et al. 
(1994), which show the apparent magnitude distribution of 
Galactic bulge red clump giants at the same low-extinction 
line of sight directions (±5°, —4°) in the OGLE microlensing 
surveys. While it is premature to use the simple analytical 
models here to draw any conclusion on the bar angle, they 
suggest that an integration of the COBE map with the star 
count data of bulge red clump giants (Stanek et al. 1997) 
can lift the degeneracy of the bar models. 



§ i/q is not the peak density ^(0) = (l — ^1 vq (cf. eq. |sj) , which 
roughly doubles by going from a a = 22.5° model to a a = 63.5° 
model. 



The non-uniqueness in deprojecting the COBE/DIRBE 
maps originates from the simple fact that a round bulge can 
mimic an end-on bar in the line-of-sight integrated light dis- 
tribution. This kind of non-uniqueness is almost lifted when 
bar-like reflection symmetries are imposed, except that the 
angle of these symmetry planes with our line of sight is still 
a free parameter, which can only be loosely constrained by 
the positivity of the deprojected density. The main results 
are summarized as follows. (1) Triaxial Galactic bar models, 
which project to identical intensity map in the absence of 
extinction, form a ID sequence as function of the bar's ori- 
entation, in agreement with numerical results of Binney & 
Gerhard (1996). (2) It is ill-conditioned to speak of any best- 
fitting bar model because the limit on the bar parameters 
comes from plausibility of the density (positivity, smooth- 
ness, and regular-looking shape and radial profile, etc.), not 
from the goodness of fit. (3) The theory predicts that COBE 
models likely follow a sequence with identical v dd, which is 
uniquely constrained by the perspective effects of a triaxial 
bar. (4) The non- uniqueness in the even part of the den- 
sity v eV en is of nature entirely different from those induced 
by shearing/stretching of an external ellipsoid (Stark et al. 
1977) or adding a konus-like structure (Gerhard & Binney 
1996). (5) Various methods to fully break the degeneracy of 
COBE bars are compared. In agreement with findings for 
extragalactic axisymmetrical systems (van den Bosch 1997, 
Romanowsky & Kochanek 1997), it is unlikely to put a very 
tight limit on the bar's parameter space with general dynam- 
ical constraints such as self-consistency and stability, but it 
is still hopeful to lift some degeneracy by fitting kinematic 
data of gas and stars with sophisticated dynamical models. 
With the increasing number of microlensing events towards 
the bulge, microlensing optical depth maps also have the 
potential of lifting the degeneracy of the models. But the 
cleanest approach to remove the ambiguity of the bar an- 
gle is to compare with star count data from the current 
bulge microlensing and variable star surveys, in particular 
the large data set on the bulge red clump giants, which pro- 
vide the extra constraints on the "depth" of the bar (Stanek 
et al. 1997). Extinction is the main limiting factor in analysis 
of both star counts and COBE maps. Finally the mass-to- 
light ratio, the last parameter of a complete mass model of 
the bar, can be determined from fitting gas/stellar kinemat- 
ics or the microlensing optical depth (Bissantz, Englmaier, 
Binney, & Gerhard 1997). 

I thank Frank van den Bosch for a helpful discussion 
in the early stage of this work, Tim de Zeeuw and Ortwin 
Gerhard for a careful reading of an early draft. 
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